Astronomy & Astrophysics manuscript no. nschneider 


©ESO 2010 


December 24, 2010 





The link between molecular cloud structure and turbulence 

N. Schneider^, S. Bontemps^, R. Simon^ V. Ossenkopf^ C. Federrath'*'^'^, R.S. Klessen^-'^, F. Motte\ Ph. Andrei J. 

Stutzki^, and C. Brunt^ 

' IRFU/SAp CEA/DSM, Laboratoire AIM CNRS - Universite Paris Diderot, 91191 Gif-sur-Yvette, France 

2 OASU/LAB-UMR5804, CNRS, Universite Bordeaux 1, 33270 Floirac, France 

^ I. Physikalisches Institut, Universitat zu Koln, Zulpicher StraBe 77, 50937 Koln, Germany 

■* Zentrum fiir Astronomic der Universitat Heidelberg, Inst, fiir Theor Astrophysik, Albert-Ueberle Str 2, 69120 Heidelberg, Germany 

^ Max-Planck-Institut fiir Astronomic, Konigstuhl 17, 691 17 Heidelberg, Germany 

* Ecole Normale Superieure de Lyon, CRAL, 69364 Lyon Cedex 07, France 

' Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Menlo Park, CA 94025, U.S.A. 
School of Physics, University of Exeter, Exeter, EX4QL, UK 

December 24, 2010 

ABSTRACT 

Aims. We aim to better understand how the spatial structure of molecular clouds is governed by turbulence. For that, we study the 
large-scale spatial distribution of low density molecular gas and search for characteristic length scales. 

Methods. We employ a 35 square degrees '""CO 1^0 molecular line survey of Cygnus X and visual extinction (Ay) maps of 17 
Galactic clouds to analyse the spatial structure using the A- variance method. This sample contains a large variety of different molecular 
cloud types with different star forming activity. 

Results. The A-variance spectra obtained from the Ay maps show differences between low-mass star-forming (SF) clouds and massive 
giant molecular clouds (GMC) in terms of shape of the spectrum and its power-law exponent p. Low-mass SF clouds have a double- 
peak structure with characteristic size scales around 1 pc (though with a large scatter around this value) and 4 pc. GMCs show no 
characteristic scale in the Ay-maps, which can partly be ascribed to a distance effect due to a larger line-of-sight (LOS) confusion. The 
A-variance for Cygnus, determined from the "CO survey, shows characteristic scales at 4 pc and 40 pc, either reflecting the filament 
structure and large-scale turbulence forcing or - for the 4 pc scale - the scale below which the '^CO 1— >0 line becomes optically 
thick. Though there are different processes that can introduce characteristic scales, i.e. geometry, decaying turbulence the transition 
scale from supersonic to subsonic turbulence (the sonic scale), line-of-sight effects and energy injection due to expanding supernova 
shells, outflows, Hll-regions, and the relative contribution of these effects strongly varies from cloud to cloud, it is remarkable that 
the resulting turbulent structure of molecular clouds shows similar characteristics. 

Key words, interstellar medium: clouds - individual objects: Cygnus X - molecules - kinematics and dynamics - Radio lines: ISM 



[ 1. Introduction 

. Turbulence 

' Over the last years, much progress has been made to theoret- 
, ically explain the observational claim that turbulence governs 
. the different phases of the Interstellar Medium, ISM (MacLow 
■ & Kle ssen l2U04l Elmegreen & Scalo l2004l Scalo & Elmegreen 
, I2004I I. Molecular line surveys of molecular clouds as a whole 
led to the detection of the linewidth-size cr oc L'Q and mass-size 
[M oc U relations (La rson [T98n i with e ~ 0.2...0. 5 and y ~ 2 
! (e.g. G oldsmith et al. [19871 He yer & Brunt l2004l Heyer et al. 
' 120091 Lombardi, Alves & Lada I2010I I. and a mass distribution 
for molecular clumps of the form of dN/dM ~M"" with a -1.5- 
1.9 (e.g. Stutzki & Gusten lT990l Kramer et al. lT998l l. The spa- 
tial structure of the emission has been characterised in terms of 
power spectra (Scalo [T9871 Stutzki et aL 119981 ). 

The scaling relations and power laws can be explained in 
the context of turbulence, where energy is injected at very large 
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' In the following, we use a for the clump mass spectral index, /? for 
the slope of the A-variance, y for the power-law exponent of the mass- 
size relation, and e for the power-law exponent of the linewidth-size 
relation. 



scales and cascades down to the smallest scales, creating eddies 
and stirring up the cloud. In particular the linewidth-size rela- 
tion of giant molecular clouds can only be reproduced if ISM 
turbulence includes sources of driving that act on large scales, 
but small-scale driving by outflows can be important in localised 
(sub-parsec) regions of larger clouds. The driving sources for 
turbulence on the largest scales can be diverse, ranging from the 
accretion of gas of extragalactic origin (Klessen & Hennebelle 
1201 Oi l to the occurance of convergent flows of atomic gas trig- 
gered by spiral density waves (Walder & Folini [19981 Vazquez- 
Semadeni et al. 120071 Heits ch et al. 120081 Hennebelle et al. 
20()8] Banerjee et al. I2009I I. supernova explosions (MacLow 
& Klesse n [20041 Dib et al. [20U9l l. or expand ing H II regions 
( Matzner 120021 Krumholz, M atzner, Mc Kee 120061 P eters e t 
al. [20081 Krumholz & Matzner [2009l Gritschneder et aL [2009l l. 
Here it is the very process of cloud formation that drives the in- 
ternal turbulence. Observations by Brunt et al. ( I2009I I indicate 
that the size-velocity dispersion of giant molecular clouds can 
only be reproduced if ISM turbulence is driven on large scales. 
Some models have investigated molecular cloud turbulence that 
is driven on small scales by internal sources such as stellar winds 
and o utflow s (Li & Nakamura [2006l Nakamura & Li [2008l Wang 
et al. 1201 Oi l but it is unlikely that these have a significant effect 
on the largest scales within clouds (Mac Low & Klessen 120041 
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Banerjee et al. 120071 Brunt et al. 12009b . All the proposed driving 
sources above will excite a considerable amount of compressive 
modes (Federrath, Klessen, Schmidt 120081, which are essential 
to trigger star formation. 

In its simplest form, the Kolmogorov-turbulence 
(Kolmogorov ll941l ) describes an incompressible, non-magnetic 
fluid in which the statistical properties of the flow are inde- 
pendent of the scale within the inertial range. This leads to 
an index of e=0.33 and to a self-similar, i.e. fractal, spatial 
structure. Other turbulence fluid models (Vazquez-Semadeni et 
al. 119971 Baflesteros-Paredes & MacLow 2002. Federrath et al. 
I2010I I. describing a supersonic, isothermal turbulent flow obtain 
e=0.5 and y=2. Hennebel le & A udit ( I2007I I and Hennebelle, 
Audit, Miville-Deschenes (12007b take a diff'erent approach and 
describe turbulence by the individual motion of cloudlets. Their 
hydrodynamic, non-magnetic, non-self-gravitating simulations 
of a turbulent 2D-phase of atomic HI gas led to a mass-size 
relation of M ~L'^ and a clump mass spectral index a ^1.7, a 
value similar to that found for CO clumps. 

These few examples show already the effort that is made 
to take into account the complexity of the ISM and to com- 
pare observations with theory. However, each model represents 
a particular physical situation with different focus on the various 
physical effects that play a role (driving source of turbulence, 
magnetic fields, self-gravity, thermal balance, etc.). In particular 
the importance of self-gravity was pointed out by, e.g. Klessen 
( 12000b who discussed its effect on various statistical characteris- 
tics, and Vazquez-Semadeni et al. (12008b who emphasized that 
a part of the observed velocity dispersion in molecular clouds is 
due to clump-scale inward motions as a consequence of gravita- 
tional collapse and not turbulence. Recently, Ballesteros-Paredes 
et al. ( 12010b have suggested that the size-linewidth-surface- 
density relation observed by Heyer et al. (12009b can be quan- 
titatively explained if the internal turbulent motions are gener- 
ated by hierarchical contraction. Their model successfully ex- 
plains the apparent virialisation of the clouds and potentially 
avoids the need for carefully balanced 'replenishment' of the 
turbulence by 'internal' driving sources. Early explanations of 
linewidths in molecular clouds (Goldreich & Kwan 11974b as 
arising from gravitational collapse were opposed (Zuckerman 
& Palmer 11974b on the basis of the untenably high star for- 
mation rates that would ensue. Ballesteros-Paredes et al. (12010b 
propose that this problem may be avoided if the higher density 
regions collapse first (by virtue of their shorter free-fall times) 
and induce star formation feedback which interrupts the global 
collapse and quickly destroys the cloud. Recent hydrodynamic 
models (Kr umhol z l2006l Krumholz et al. l2007l Offner et al. 
120091 Bate 12009b stress the importance of radiative feedback 
processes from stars that can drastically reduce the total number 
of stars formed, although it has been found that fragmentation is 
never completely suppressed (Peters et al. 120 1 Oal 120 1 Obi |20 1 Ocl ) 
Similar holds for the influence of magnetic fields in massive star- 
forming regions (Hennebelle et al. 120101 Peters et al. l2010d] l. 

These 'deviations' from the idealized scenario of a purely 
turbulent medium may cause changes of the scaling laws. This 
could explain, for example, why many molecular line studies 
failed to detect a linewidth-size relation in molecular clumps 
(e.g. Loren [19891 Simon et al. [20011 Schneider & Brooks 
12004b . In any case, it is difficult to transform model results 
into 'observables' such as the 3D- or projected (on the plane 
of the sky) emission and velocity structure of an observed 
molecular line. Ballesteros-Paredes & Mac Low (2002) argued 
that the apparent mass-size relation is due to the observations 
and/or the analysis process: any cutoff in column density, due 



to the limited dynamic range and/or the minimum intensity 
level defined by clump finding algorithms automatically implies 
M cc Lr. Even more problematic is the velocity superposition 
of many clumps in the line-of-sight when a low-density tracer 
like '"CO is used. High-density tracers like CS or N2H^ 
provide a better clump seperation but imply a bias in analysing 
the full cloud structure (only the densest structures are identi- 
fied). A survey using the '^CO line may be the best compromise. 

Molecular cloud structure 

The observed spatial structure of molecular clouds depends on 
the tracer used to image the clouds. Low-J rotational transitions 
of '-CO or visual extincion maps are sensitive to low densities 
and show a very complex structure that - with a change to 
higher angular resolution - breaks up into substructures, which 
appear similar on almost all scales (e.g., Falgarone et al. 119911 
Stutzki [2001b . The concept of 'self-similarity' as a geometrical 
property was introduced (Elmegreen & Falgarone |1996l ) and this 
gas component is well described as a fractal-Brownian-motion 
(fBM) structure (Stutzki et al. 119981 ). In contrast, gas at much 
higher column density, best traced by dust (sub)mm-continuum 
emission or molecular line high density tracers like CS, shows 
no fractal structure (Falgarone et al. 120041) . In addition, self- 
similarity may break down on the smallest scales (<0.1 pc), 
relevant for star formation (e.g. Goodman et al. 119981 Dib et al. 
[20081 Federrath et al. [20T0l l. 

Describing molecular cloud structure 

Statistical tools have been developed to analyse the observed 
spatial and kinematic cloud structure, which turns out to be 
extremely complex. Amongst the most commonly used tools to 
describe structure are probability distribution functions of line 
intensities (Padoan et al. 11997b . the A-variance analysis of inten- 
sities and velocities (Stutzki et al. 119981 Ossenkopf 12002b . the 
principal component analysis (Heyer & Schloerb 119971 Brunt 
& Heyer I2002al I2002bl ). and structure functions (Federrath et 
al. ZOrOj). In the current study, we use the A-variance method, 
applied to a large '''CO 1^0 survey obtained with the FCRAC0 
14m telescope of the Cygnus X region and on a data set of 
17 near-IR extinction maps of molecular clouds derived from 
2MASS. This analysis will help to further develop and refine 
models and cloud structure characterisation techniques. We are 
particularly interested to see whether there is a dependence 
of cloud characteristics (size scales, value of the power- 
law slope j6, form of the A-variance spectrum) and cloud type 
(low- and high-mass star-forming molecular clouds, diffuse gas). 

Cygnus X 

The Cygnus X region is the richest and most massive com- 
plex of high-mass star formation at a distance lower than 3 kpc. 
Schneider et al. ( 120061 12007b showed that the Cygnus X region 
constitutes a large scale network of GMCs at a common distance 
of ~1.7 kpc. Signposts of recent and ongoing (high-mass) star 
formation (H ii regions, outflow activity, IR-sources such as S 106 
IR, D R21, W 75N, and GL2591) are ubiquitous. See Schneider 
et al. ( 12006b or Reipurth & Schneider (12008b for a review on 
Cygnus X. The region also contains several OB clusters (Cyg 
OB 1-4, 6, and 8) including the richest known OB cluster of 
the Galaxy, Cyg OB2 (Knodlseder [2000l Uyaniker et al llUOTT l. 
Their energy injection from stellar winds and radiation could be 
a driving source for turbulence. 
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Fig. 1. Line integrated map of Cygnus X in '^CO 1^0 emission with the intensity (-1 to 16 K km s ') coded in color. Crosses 
indicate the location of thermal HII regions, the well-known star-forming regions DR21 and SI 06 are labelled. 



2. Observations 

2.1. FCRAO '^CO 1^0 survey 

The data presented here were obtained using the FCRAO partly 
in remote operation from late 2003 to early 2006. The sum- 
mer period was excluded for observations due to weather con- 
straints. The molecular lines of '^CO 1-^0 at 110.201 GHz and 
C"^0 1^0 at 109.782 GHz were observed simultaneously. The 
beamwidth of the FCRAO at 110 GHz is 46". 

A total of 35 square degrees in Cygnus X was covered 
in the '^CO and C'^O l->0 lines at 46" angular resolution 
(Simon et at, in prep.). This survey is significantly larger and at 
higher angular resolution than the '''CO 2— >1 map obtained with 
KOSMAB presented by Schneider et al. (2006) . It has a spec- 
tral dynamic range (maximum observed linewidth over velocity 
resolution) of ~70 and a spatial dynamic range (map size over 
angular resolution) of ~230. This is much larger than what is 
usually obtained for molecular cloud maps (Kramer et al. |1998t . 
The highest column density regions traced in the FCRAO '"'CO 
1— >0 survey were then observed in the CS 2-^1 and N2H^ 1^0 
lines (~15 square degrees) using the FCRAO. 

Pointing and calibration were checked regularly at the start 
of the Cygnus observing time interval and after transit of Cygnus 
(no observations were performed at elevations higher than 75°). 

' Cologne Observatory for Submm- Astronomy 



Pointing sources were SiO masers of evolved stars, i.e., ;t'-Cyg, 
R-Leo and T-Cep, depending on observing time. The calibration 
was checked regularly on the position of peak emission in DR21 
and found to be consistent within 10%. 

We used the single sideband focal plane array receiver 
SEQUOIA (Second Quabbin Optical Imaging Array; Erickson 
et al. 1 19991 1 designed for the 85-115.6 GHz range, yielding a 
mean receiver noise temperature of 60 K. SEQUOIA is a 4x4 
dual-polarization array with a separation between elements of 
88" on the sky. The receiver was used in combination with a 
dual channel correlator (DCC), configured to a bandwidth of 25 
MHz, using 1024 channels, corresponding to a velocity sampling 
of 0.066 km s'. 

We used the on-the-fly mapping mode (OTF) in which each 
map consisted of a block of 20'xl0' (30 minutes integration 
time including overheads). The data quality is enhanced because 
pixel-to-pixel variations are averaged in the final regridded data 
set. Since at the map edges the coverage is less dense, and thus 
leads to an increase in noise, adjacent maps were overlapping to 
cover the region at uniform density. 

In total, ~800,000 spectra were produced on a 22". 5 sam- 
pled grid covering an area of 35 deg-. The data have a mean 1 
CTrms rms noisc level of ~0.2 K per 0.06 km s ' channel on a T^ 
antenna temperature scale, i.e., not corrected for the main beam 
efficiency of ~0.48. More details on the observing procedure, 
regridding and the rms-map of '^CO are found in Appendix A. 
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2.2. Near-Infrared extinction maps 

The basic idea of deriving extinction maps from near-infrared 
(near-IR) data is that the average reddening of background stars, 
as observed in the near-IR bands (usually J, H and K), is used to 
estimate the total column density of dust along the line-of-sight 
and thus the total Ay in that direction (e.g. Lada et al. 119941 ). 

We developed a new technique (AvMAP) to produce in an 
efficient way extinction maps from the 2MASS0 catalog fol- 
lowing the schemes described in Lombardi et al. C200T| and 
Cambresy et al. (I2002I I. but with the improvements that (i) it 
is a fast algorithm which can easily handle large fields with 
millions of sources, (ii) include a calibration of the extinction 
which does not require an off-field, and (iii) it is a method to 
estimate the density of foreground stars which follows the spa- 
tial distribution of stars in the Galaxy by using the predictions 
of the stellar population model of Besangon (Robin et al. 120031 
see http://model.obs-besancon.fr/) to filter out foreground stars. 
Thus, Ay-maps can be used up to a distance of 3 kpc at an an- 
gular resolution of typically 2'. A detailed description of the 
method is given in the appendix B. 

3. The A-Variance as a method of structure analysis 

3.1. Definition of A-variance 

The A-variance is a method to analyse the structure and correla- 
tion in multi-dimensional scalar functions. It was first introduced 
by Stutzki et al. ( II998I ). who expanded the 1-D Allan variance 
concept (Allan 1 196 61) to two and more dimensions in order to 
characterise the observed 2-D structure in molecular cloud im- 
ages. 

The A-varianc^ measures the amount of structure on a given 
scale L in a map s (2-dimensional scalar function) by filtering 
this map with a spherically symmetric wavelet Ql. 

The A-variance thus probes the variation of the intensity s 
over a length L (called lag), i.e., the amount of structural varia- 
tion on that scale. The filter function can vary between a smooth 
or a step-shaped filter function ('Mexican hat' or 'French hat'). 
Ossenkopf, Krips, Stutzki (I2008al ) analysed systematically (us- 
ing simulations) the influence of the shape of the filter function 
and arrived to the conclusion that the Mexican hat filter, with 
an annulus-to-core-diameter ratio of about 1.5, provided the best 
results for a clear detection of pronounced scales. 

The observed structure can be well mimicked by a fractional 
Brownian motion (fBm) structure (which has a power-law spec- 
trum) in the context of fractal images. For any 2-D image with a 
power spectrum P{k) oc \k\^^, in which k is the spatial frequency, 
the 2-D A-variance varies as cr^ oc lP^^ for < y6 < 6 (see 
Stutzki et al. 1998 for more details). The range of L is given by 
the upper and lower spatial frequency cutoff's k^ ~ (Ajc)"' and 
ki » (nAxy^ (Ax is the sampling grid of a map with pixels) : 
(Inku)-^ <L< (2nk,)-\ 



The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 

^ All A-variance calculations were performed using the IDL- 
based routine deltavar provided by V. Ossenkopf and available at 
www.astro.uni-koeln.de/~ossk/ftpspace/deltavar 
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Fig. 2. The A-variance spectrum of the line integrated map of the 
entire complex in the '^CO emission. The slope /? is determined 
from a linear fit for data between lags 1' and 8' (indicated by the 
grey shaded area). The dotted line gives the beamsize (45")- 

The advantage of the A-variance compared to the power- 
spectrum is that Fourier-based methods, like the determintation 
of the power-spectrum, can create artifical structures by edge 
and gridding effects. This is not the case for the A-variance, 
the size of the image only limits the largest lag probed and the 
shortest lag will be close to the angular resolution. In addition, 
noise contributions in the map can be separated from the intrin- 
sic structure of the cloud because purely white noise has fi-Q 
(flicker noise would lead to /?=1). Weighting the image with the 
inverse noise function (l/cr, ,,, ,) enables us to distinguish the vari- 
able noise from real small scale structure (see Bensch et al. 2001 
for more details). At the smallest scales, the A-variance spectrum 
is dominated by the blurring from the finite telescope beam and 
radiometric noise. 

4. Results 

4.1. FCRAO^^CO 1^0 data 

Figure [U shows the line intergrated map of Cygnus X in '^CO 
1— >0 emission, tracing gas densities around 3x10"^ cm""*. The 
main features seen in Cygnus X are high-contrast filaments, 
mainly in the northern part of the map (DR21), and more dif- 
fuse emission with some higher density clouds in the south. 
Schneider et al. (120061 ) named these two regions 'Cygnus North' 
and 'Cygnus South'. The molecular line emission from those 
clouds related to the Cygnus X cloud complex covers a velocity 
range from roughly -15 to +20 km s ' (Schneider et al. |2006]). 
However, some emission at velocities larger than about 4 km s ' 
(Piepenbrink & Wendker [T988l Schneider et al. 12007] ) is due to 
the 'Great Cygnus Rift', a low-density gas feature at a distance 
of around 600 pc. It is responsible for a visual extinction of a 
few Ay (Dickel [19691 Schneider et al. 12007] ). In the figure, sev- 
eral arc or shell-hke structures are visible, in particular in the 
northern part of the map where a cavity devoid of molecular gas 
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(center at offsets 10', 130') is surrounded by clouds that may 
have formed as swept-up material from the HII regions and the 
influence of the OB2 cluster In Appendix C, we show two ad- 
ditional figures, FCRAO '^CO 1^0 data overlaid on a map of 
radio continuum emission at 1420 MHz from the CGPS0 and on 
mid-IR emission from MSXQ in order to better depict the distri- 
bution of neutral/ionized and warm/cold gas. 

In order to show the variation in the amount of structure as 
a function of length (lag) in the clouds, we apply the A-variance 
method to the line integrated '^'CO 1— >0 map. This is the first 
time that a map of ~35 square degrees is analysed with the A- 
variance. Typical map sizes in previous studies (Bensch et al. 
120011 Sun et al. 12006b were in the range of 1 to 10 square de- 
grees. We profited from the (1/rms) weighting feature imple- 
mented in the A-variance algorithm developed by Ossenkopf, 
Krips, Stutzki (2008a), so that areas with higher noise (in our 
case the map edges) contribute less to the A-variance than ar- 
eas with low, uniform noise. This improves the accuracy of the 
results with respect to distinguishing real signal from noise. 

Figure |2] shows the A-variance spectrum for the '^CO 1— >0 
map. The error of each cr^(i) is determined from the Poisson 
statistics. In Appendix D, we explain in more detail the er- 
ror calculation. The plot shows that with increasing lag L, the 
A-variance first increases and reaches its peak value at Lpeak, 
i.e. 8' corresponding to 4 pc for a distance of 1.7 kpc. With 
increasing lag, a second peak is observed at 80' (40 pc) be- 
fore the A-variance decreases approximately with L"^. Because 
this spectrurr0 exhibits well defined characteristic scales, it does 
not represent a fully self-similar structure which would be the 
signpost of fully developed supersonic turbulence (Mac Low & 
Ossenkopf I2000I I and which would produce a pure power law. 

The existence of characteristic scales in the A-variance spec- 
trum could be due to (i) the true presence of preferred geometric 
scales in the system (length and width of filaments for example), 
(ii) decaying turbulence (Mac Low & Ossenkopf 12000b . (iii) en- 
ergy injection from external or internal sources (larger-scale ac- 
cretion, spiral density waves, supernovae, outflows and winds, 
expanding HII regions, etc.), or by (iv) superposition of differ- 
ent clouds along the line-of-sight (note, however, that this effects 
does not create characteristic scales but simply leads to several 
peaks in the A-variance spectrum). It is therefore not straight- 
forward to interpret the A-variance spectrum of Cygnus since 
several of these effects may play a role. 

We can, however, argue that the characteristic scales seen in 
the A-variance spectrum are not caused by spurious line of sight 
effects. It is known that at least two groups of molecular clouds 
at distances of 600 pc ('Great Rift') and 1.7 kpc (see above) are 
found in the direction to Cygnus X. But it was shown and dis- 
cussed in Schneider et al. (12007b that the column density of the 
Rift gas is very low (Ay <3'") and not well traced by '^CO emis- 
sion. It is thus well justified to assume that the whole A-variance 
spectrum would be due to emission features arising from the 
clouds in 1.7 kpc distance. 

Characteristic scales due to energy injection are plausi- 
ble since this region is strongly influenced by radiation from 
the Cygnus OB2/OB9 clusters and a large number of outflow 
sources. The typical sizes of HII regions in the field (see fig- 
ure C.l in Appendix C) are of order of a few pc and larger. 



* Canadian Galactic Plane Survey 
^ Midcourse Space Experiment 

^ None of the scales corresponds to a structure introduced by the 
observational procedure because individual OTF-maps have a size of 
10'x20'. 



thus fitting to explain the origin of the observed scale of 4pc. 
Outflows, however, act on much smaller scales. Nakamura & Li 
(2007) report a characteristic scale of 0.3 pc from their simula- 
tions of outflow driven turbulence. This appears consistent with 
the scales observed in low-mass star forming regions (e.g. as 
seen in velocity channel maps in NGC 1333, Quillen et al. 120051 
Curtis et al. 12010b . However, Cygnus X forms a large number of 
high-mass stars which are expected to drive more energetic and 
most likely larger outflows compared to their model. A conclu- 
sive answer, therefore, requires more detailed numerical studies 
with focus on high-mass star forming regions. We also note that 
only a few signatures of SN-shells are found in Cygnus X and 
we conclude that the rate of SN-explosions in this area is very 
low (see overview in Reipurth & Schneider 120081 ). 

Another possible explanation for the features in the observed 
spectrum is that with '''CO we indeed trace signatures of true 
morphological characteristics of the cloud. In concert with the 
visual inspection of the region (Figure[T]i it is appealing to iden- 
tify the first peak around 4 pc with small-scale filaments as also 
seen in many other nearby clouds (see Section |4. 2. 2b . The sec- 
ond peak corresponds to the large scale filamentary structure that 
dominates the global appearance Cygnus X. 

However, there is another caveat to notice. An alternative 
explanation to the 4 pc scale could be that it is due to radia- 
tive transfer effects in turbulent media (Padoan et al. 120001 ). As 
'''CO 1— »0 turns optically thick at column densities above about 
N(H2)=10^^ cm~2 (Ossenkopf l2002b . this tracer is not capable of 
resolving any denser structure and all peaks are bluiTed out. The 
'^CO observations may thus only provide a good reflection of 
the structure measured up to about Ay =5"'. The peak in the '■'CO 
A-variance basically shows the typcial size of the regions where 
the tracer turns optically thick, i.e. the '^CO "photosphere" of the 
dense regions. The exact location of the peaks on the A-variance 
spectrum gives some measure for the turbulence in the molecular 
cloud. They indicate a transition from more coherent structures, 
where the '^CO emission falls into a narrow velocity interval, to 
a wide-spread turbulent distribution leading to a dilution of the 
optical depth in very broad lines. The details of this transition 
and its effects on the physical interpretation of the spectrum will 
be studied in a subsequent paper based on the A-variance anal- 
ysis of the velocity structure, i.e. studying individual velocity 
channel maps. 

The power spectral index /3 is typically determined for lags 
between the beamsize (indicated by a dashed line in Fig.|2| and 
the first significant structure component (8'). The value /3-2.53 is 
at the lower end of typical values for molecular clouds (see, e.g., 
Bensch et al. ( I200II ) for a complilation of clouds with indices be- 
tween 2.5 and 3.3). This implies more structure on smaller scales 
and thus that the lower column density gas has more structure. 
This was already pointed out by Falgarone et al. (120041) who 
studied the diffuse component of the interstellar medium, find- 
ing highly dynamic and fractal structures. 

4.2. Extinction maps 

4.2.1. A-variance from extinction maps 

Figures |3] to |6] show extinction maps of Galactic molecular 
clouds, obtained with the method described in Appendix B, and 
the corresponding A-variance. For regions where '''CO data was 
available from the literature, we aditionally plotted the obtained 
A-variance spectra. We determined the A-variance without rms- 
weighting in the IDL routine 'deltavar', using the Mexican-hat 
filter with a diameter ratio of 1.5, recommended by Ossenkopf, 
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Fig. 3. Extinction maps of Galactic molecular clouds (left) and corresponding A- variance spectrum with errors (right). See Appendix 
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Krips and Stutzki (!2008ab . All maps are extended and include all 
major emission regions (note that all A-variance spectra cover 
lags well below the map extent). However, "empty" regions also 
contribute to the A-variance appearing as large coherent struc- 
tures so that they produce a decay of cr^ at large lags which is 



shallower than the L ^ decay that would represent a zero corre- 
lation at those scales (Ossenkopf, Krips, Stutzki l2i308al l. 

The regions that were analysed comprise close-by (distance 
<1 kpc), low-mass star-forming regions like Perseus, Taurus or 
the Pipe Nebula, as well as more distant (up to 3 kpc) high-mass 
star-forming regions like Cygnus, Rosette or W3. 
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A number of spectra show a simple rising curve (e.g. 
NGC2264, Coalsack, Vela), some including a noise de- 
crease plus a contribution due to empty regions (W3, Taurus, 
Chameleon). The majority of spectra, however, have more com- 
plex shapes, (e.g. Lupus, MonR2, MonOBl, IC 5146, Pipe) 
with positive or negative curvature and/or two more or less 
pronounced peaks. Interestingly, the most filamentary struc- 
tured clouds, partly with very thin, elongated features (Taurus, 



Chameleon) have the most complex double-peak curves and 
show low values of /3. 

It is tempting to interprete, in a purely geometrical view, the 
two peaks in the A-variance spectrum (<1 pc and ~4 pc) as the 
width and length of filaments. However, recent studies of the 
filamentary structure in the Aquila and Polaris regions with the 
Herschel satellite (Andre et al. 120 101 Men'shchikov et al. 120 1 Oi l, 
obtained a much smaller typical width of the filaments. A de- 
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Fig. 5. Extinction maps of Galactic molecular clouds (left) and corresponding A-variance (right). The red curve for Cygnus repre- 
sents the '^CO 1^0 data from the FCRAO survey. See Fig. 3 for further explanations. 



tailed analysis of the radial profiles of the filaments observed 
with Herschel in Polaris, Aquila, and ICS 146 resulted in a char- 
acteristic filament width of about 0. 1 pc was (Arzoumanian et 
al. (in prep.). The shape of the A-variance and the value of /3 
weakly depend on the size of the selected region. Subregions in 
a map (e.g. for Vela, where only the high column density re- 
gion was selected from the Ay map), generally produce slightly 
steeper spectra with higher values of /3. Large maps showing ex- 
tended filamentary structures produce flatter spectra (e.g. Taurus, 
Chameleon, Pipe). 



4.2.2. Comparison to '^CO data 

For regions where we had literature values or own data, we in- 
cluded the A-variance determined from '^CO, i.e. Perseus, NGC 
2264, Orion A+B, Rosette, Cygnus, and Mon R2. The spectra 
of the massive GMCs Cygnus and Rosette look significantly dif- 
ferent compared to the spectra derived from the extinction maps. 
For Cygnus, the curve determined from the Ay-maps represents 
a fully self-similar scale distribution, covering all scales (since 
it is a very large and extended map), while the '^CO spectrum 
shows the double-peak feature already discussed in Sec. 14. II We 
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attribute the two peaks from the '^CO curve at ~4 and ~40 pc to 
the typical filament-sizes in Cygnus X, arising from clouds at a 
distance of 1.7 kpc and/or optical depths effects of the '^'CO line 
emission. However, these two characteristic scales are not appar- 
ent in the A-variance spectrum of the extinction map. A possi- 
ble explanation is that since the Ay map comprises all emission 
features along the line-of-sight, including the one from the low- 
density Cygnus Rift in 600 pc distance, the resulting spectrum 
appears indeed self-similar since it is a composite of diffuse, 
fractal gas and denser structures. It is thus possible that by using 



the '^CO line, low density larger scale structures that do not emit 
in '^CO around ~10'-~100' are partly filtered out. However, the 
two spectra show a similar increase of the A-variance (up to a 
scale of about 10') with the slope of the extinction map being 
slightly shallower. 

The Rosette spectrum has the same tendency to show a 
peaked structure in '''CO with a plateau around 8 pc and then 
a slight decrease, while the extinction map has no characteristic 
scale. However, for small lags, the spectra show a very differ- 
ent behaviour Since the '^CO spectrum is not affected by noise 
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for small lags, the spectrum increases in a similar way like for 
Cygnus (up to 20')- The spectrum from the Ay-map, on the other 
hand, is limited by noise for small lags and is thus not directly 
comparable to the '^CO spectrum. The opposite is true for large 
lags. The 8 pc scale may again indicate the scale when the '^'CO 
line becomes optically thick. The larger scale for Rosette com- 
pared to Cygnus (4 pc) could be due to the fact that the average 
density of clumps in Cygnus (Schneider et al. 120061 ) is higher 
than in Rosette (Williams et al. 119951 ) and thus the ^^CO line is 
saturated on a smaller scale. 

For the other sources (Perseus, Orion A/B, Mon R2, and 
NGC2264), we obtained the A-variance from Fig. 10 in Bensch 
et al. (l200TT l. based on '^CO l->0 data from the Bell Labs 7m 
telescope survey (Bally et al., priv. communication). The form 
of the spectrum is rather similar to the ones obtained from the 
Ay-maps but the location of peaks and the values of fi slightly 
differ. It seems that the angular resolution of the maps has an 
influence on the shape of the A-variance since its form for the 
extinction maps (at 2' resolution) correspond well to the 1.7' 
Bell Labs data but differ from that of the 50" resolution FCRAO 
surveys. However, Padoan et al. (12003b found for Taurus and 
Perseus that the structure function of '^CO follows a power-law 
for linear scales between 0.3 and 3 pc, similar to our finding 
that a first characterestic scale is seen around 2.5 pc. A more re- 
cent study of Brunt (1201 Oi l of Taurus, comparing the power spec- 
tra determined from '^CO and Ay showed that they are almost 
identical and that there is a break in the column density power 
spectrum around 1 pc. Below a wavenumber corresponding to 
a wavelength of ~ 1 pc. Brunt ( 120101 ) found a power spectrum 
slope of 2.1 (similar to the delta variance slope found here). At 
wavenumbers above that (smaller spatial scales), a steeper spec- 
trum was seen, of index 3.1. This power spectrum break is as- 
sociated with anisotropy in the column density structure caused 
by repeated filaments (Hartmann 120021 ). possibly generated by 
gravitational collapse along magnetic field lines. 

4.2.3. The A-variance for all clouds 

Figure [T] shows the A-variance spectra for all clouds in this 
study obtained from the Ay-maps. At small scales (below 1 
pc) most sources show a non-constant spectral index steepen- 
ing towards the resolution limit. This is consistent with decay- 
ing turbulence dissipated at small scales (Ossenkopf &Mac Low 
120021) but also with driven turbulence at small lags (Federrath 
et al. 120091 ). However, there are three exceptions: Chameleon 
and Taurus show an intermediate peak at 0.3-0.4 pc and Perseus 
shows no clear indication of a steepening. The intermediate peak 
might mean that the extinction map is affected by a separate, 
more distant, component that is actually dominated by struc- 
tures larger than assigned in the plot (see also discussion be- 
low). Alternatively, it could be produced by a systematic struc- 
ture of the detected size that affects the turbulence in the cloud. 
Candidates for such structures are SN-shells. Expanding ioniza- 
tion fronts from OB associations impact the cloud structure as 
well. For example, it is known that Lupus is influenced by a sub- 
group of the Sco OB2 association (Tachihara et al. 120011 ) both by 
past SN explosions and present OB stars and their HII regions. 
Cygnus is exposed to the very massive OB2 cluster (Knodlseder 
et al. 120001) but lacks a significant number of SN shells. Orion 
A and B are influenced by stellar wind driven compression cen- 
tered on Ori OB lb. 

At scales above 1 pc nearly all low-mass SF clouds show 
a characteristic size scale as a peak of the A-variance spectrum 
(see Sec. 14.2.41 ). i.e. Cor. Australis, Taurus, Perseus, Chameleon, 



Pipe) show a common peak scale at 2.5-4.5pc. This indicates 
the scale of the physical process governing the structure forma- 
tion. This could e.g. be the scale at which a large-scale SN shock 
sweeping through the diffuse medium is broken at dense clouds, 
turning the systematic velocity into turbulence. The GMCs, on 
the other hand, show no break of the self-similar behaviour 
at all up to the largest scales mapped. The Rosette is com- 
pletely dominated by structure sizes close to the map size. At 
the largest. Galactic scales, energy injection due to, e.g., spiral 
density waves should be visible on a kpc size, well above our 
limits. 

4.2.4. Are there characteristic size scales depending on the 
type of cloud ? 

Although the A-variance determined from extinction maps 
shows interesting particularities for different clouds, one has to 
note that there can be large differences in the value of fi within 
each individual complex, as was shown by Sun et al. (2006b 
for Perseus. We thus have to be careful in interpretating dif- 
ferent values of /? and the shape of the spectrum depending on 
the type of cloud and tracer. In Table 1, we listed fi for differ- 
ent gas phases, from diffuse HI gas to low-and high-mass star- 
forming regions, obtained from maps with different tracers (HI, 
isotopomeric CO lines, extinction maps and dust continuum). In 
total, the values range from 2.07 (W3) to 3.12 (Vela) and are 
thus in the typical range of values determined for other stud- 
ies (e.g. Bensch et al. 12001b . However, if we exclude the most 
extreme data points for W3 and Vela and the very dispersed re- 
sults for Perseus, the range of y6-values is much more narrow, 
only 2.16 to 2.76. But in any case, no systematic trend can be 
seen between the value of fi and the type of cloud/gas phase. 
However, all /? values determined from '^CO are systematically 
higher than the ones obtained from the extinction maps. In ad- 
dition, based only on '^CO data, there is a slight tendency that 
high-mass SF regions have lower values of fi than low-mass SF 
regions. A similar trend is reported by Kainulainen et al. (120091 ) 
in the column density PDF which deviates from the log-normal 
form and develops power-law tails in massive SF regions. We do 
not see a trend with distance (since shorter distances to a cloud 
resolve smaller scale structures and thus fi should be higher). 

We also investigated the presence of a characteristic scale 
in the spectra with one or two peaks, i.e. we listed in Table 1 
the parsec-scale where a first and possible second peak were ob- 
served. Interestingly, all clouds that show a double-peak struc- 
ture are low-mass SF regions in galactic regions that are not 
(much) affected by line-of-sight crowding. The more massive 
and denser high-mass star-forming clouds show no clear double- 
peak structure. This can partly be due to the fact that most of 
them are further away (W3, Cygnus, Rosette are all at distances 
>1.5 kpc) and thus the extinction maps are sensitive for all 
gas along the line-of-sight and therefore to all structure scales. 
However, other massive clouds like NGC2264, the Orion clouds 
and Vela are closer (<1 kpc) and do not show a double-peak 
structure either. The first peak for low-mass SF clouds occurs at 
size scales that are characteristic for molcular clumps, i.e. 0.4 to 
1.5 pc (with one higher value for Perseus of 2.5 pc). The sec- 
ond peak is found between 2.9 and 4.6 pc. Blitz & Williams 
( 119971 ) found a characteristic size scale of 0.25-0.5 pc indicat- 
ing a change in structure in Taurus, conform with our finding for 
Taurus. A small sample of clouds (Orion A, p Oph, Taurus, and 
L1512) was investigated by Falgarone et al. ( 12004b using '■'CO 
1— >0 emission. They found a typical size scale of around 0.5 
pc for the width and 1 to 18 pc for the lengths of filaments. Our 
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Fig. 7. A-variance spectra of all sources where we obtained Ay maps in this study. The Lag scale is in parsec. 



sample is more homogeneous in its values and we may trace here 
two different physical processes, one leading to the formation of 
clumps (<1 pc scale) and one to the formation of filamentary 
structures (4 pc scale as the lengths of the filaments). However, 
the typical width of filaments is <1 pc as well which compli- 
cates the picture. These scales may thus be associated with tur- 
bulence driven on small scales and turbulence driven on larger 
scales. This is an intersting result since the physical properties 
and the driving vary from cloud to cloud. It might be, for ex- 
ample, that energy injection from outflows is more important for 
Perseus than for Taurus (the latter contains much fewer sources) 
and thus leads to a slightly different characterstic scale since tur- 
bulence driven by outflows acts on smaller scales than turbulence 
driven by large scale events like SN explosions. However, differ- 
ent modes of turbulence forcing can also lead to significantly dif- 
ferent A-variance spectra (Federrrath et al. |2009l l and turbulence 
statistics in the ISM (Federrath et al. 120 1 Oi l, and the same is true 
for strongly self-gravitating systems (Klessen 120001 Ossenkopf 
et al. 120011 Kainulainen et al. I2009I I. The picture might be dif- 
ferent for massive GMCs that show no characteristic scales in 
Ay maps and none in '^CO maps, if we assume that the char- 
acteristic scales around a few parsec arise from optical depths 
effects. Heyer et al. (12006b for example found no difference in 
the scaling coefficent for the velocity structure function (see, 
e.g.. Brunt l2003l for details) for the Rosette Molecular Cloud, 
an active GMC forming massive stars and strongly influenced 
by an expanding HII region, as well as G216-2.5, an even more 
massive cloud but with very little low-mass SF. They conclude 
that energy input at large scales sustains the global turbulence 
in molecular clouds while effects like expanding HII regions or 
outflows act locally in modifying the turbulent structure. 



It should also be noted that our resolution is not high enough 
to resolve molecular cores. The typical size of massive dense 
cores in high-mass star-forming regions is O.I pc (e.g. Motte et 
al. 12007b and they typically show sub-fragmentation (Bontemps 
et al. 12010b so that they are not cores in the sense of single star- 
forming units. For low-mass SF regions, like rho Ophiuchi for 
example, the typical size of a pre-stellar core is well below the 
0.1 pc scale (e.g. Motte, Andre, Neri II998l l. In any case, our 
resolution is above 0. 1 pc. It would be interesting to apply the 
A-variance on large maps at high angular resolution to resolve 
cores, clump and filaments at the same time. This will be pos- 
sible with the large-scale imaging at far-infrared wavelengths 
(70 to 500 fim) of low-and high-mass star forming regions per- 
formed with the Herschel satellite (Andre et al. 12010! Motte et 
al. l20I0i Molinari et al. 120 1 Ob . providing angular resolutions be- 
tween 6" and 40". 



4.3. A-variance and turbulence models 

To arrive at a better physical understanding of interstellar turbu- 
lence and gain deeper insight into the processes that govern the 
formation and evolution of molecular clouds it is important to 
compare the observational data with theoretical and numerical 
models of turbulent flows. This is useful, because models with 
different parameters (purely hydrodynamic or with magnetic 
field, with or without self-gravity, driven or decaying turbulence, 
different driving sources, etc.) will produce different A-variance 
spectra and thus may allow us to distinguish between differ- 
ent theories. For example, the A-variance approach was used 
by Mac Low & Ossenkopf ( 120001) to characterise the density 
and the velocity structure of interstellar turbulence simulations. 
They showed that driven, supersonic hydrodynamic turbulence 
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Table 1. Observed values of the power-law exponent /3 (column 4) in different phases of the ISM (from diffuse HI gas to GMCs). 
Values for /3 obtained from our extinction maps are indicated in bold. Column 5 and 6 indicate the size scale in parsec where spectra 
show peak values. Values in parethesis are less reliable, the spectrum shows not a prominent peak there. Values given only for peak 
2 indicate that there is only one prominent peak. Note that for Orion, we determined /3 individually for the subregions Orion A and 
B (for direct comparison to the '''CO data) and for the whole complex, which is also shown in Fig. 4. 



Region 


Distance 


Tracer 




peak 1 


peak 2 


Reference 




[pc] 






[pc] 


[pc] 




Diffuse HI 




HI 


2.6-2.8 






Green et al. (1993) 


North Polar Spur 




HI 


3.6 






Miville-Deschenes (2003) 


Diffuse Clouds 














Polaris Flare 


150 


"CO 1^0 


2.77 






Stutzkiet al. (1998) 


Polaris Flare 




"CO 2-^1 


2.76 






Falgarone et al. (1998), Stutzki et al. (1998) 


low-mass SF cloud 


Perseus 


300 


"CO 1^0 


3.07 


2.7 


(>10) 


Bensch et al.(2001) 


Perseus 




"CO 1^0 


3.09 






Sun et al.(2006) 


Perseus 




"CO 2^1 


3.03 






Sun et al.(2006) 


Perseus 




i^CO 1^0 


3.08 






Sun et al.(2006) 


Perseus 




'^CO 3-»2 


3.15 






Sun et al.(2006) 


Perseus 




extinction 


2.55 






Sun et al.(2006) 


Perseus 




extinction 


2.16 


2.5 


(>10) 




Taurus 


140 


extinction 


2.20 


0.5 


3.5 




Chameleon 


160 


extinction 


2.25 


0.4 


4.6 




Pipe Nebula 


140 


extinction 


2.36 


1.5 


3.8 




IC 5146 


400 


extinction 


2.37 


(1.4) 


(>4.5) 




Lupus 


100 


extinction 


2.49 


(1.5) 


(>10) 




Coalsack 


150 


extinction 


2.57 








Corona Australis 


170 


extinction 


2.65 




2.9 




GMC 


W3 


2200 


extinction 


2.07 




11 




Cygnus 


1700 


"CO 1^0 


2.47 


4 


40 


this study 


Cygnus 




extinction 


2.42 








Orion A 


450 


"CO 1-^0 


2.54 




6 


Bensch et al. (2001) 


Orion A 




extinction 


2.49 








Orion B 


450 


"CO 1^0 


2.68 
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(with or without magnetic field) can maintain a well-defined, 
self-similar behaviour at scales below the driving scale. Klessen 
Heitsch, Mac Low (12000b . Heitsch, Mac Low, Klessen (12001b . 
and Ossenkopf, Krips & Stutzki (2008b) showed that energy in- 
jection on a particular scale does not create density enhance- 
ments on that scale but on one which is 20% to 25% smaller. 
Finding pure power-law like A-variance spectra in molecular 
cloud data therefore indicates that the observed turbulence is 
driven from the outside. We note that virtually all theoretical 
models that applied the A-variance are based on Gaussian ran- 
dom fields or fractal Brownian motion and yield slopes that are 
well within the observed range (ai2. 3-3.1, see Table 1). It is pos- 
sible that different physical large-scale driving mechanisms (spi- 
ral density waves, supernovae, expanding HII regions) may lead 
to different spectral behaviour 

On the other hand, there are various processes that can im- 
print a characteristic scale to molecular cloud structure. For ex- 
ample, it was shown by Mac Low & Ossenkopf ( 12000b that 
freely decaying turbulence causes a change of the slope, i.e. a 



deviation from self-similarity. The only distinctive scale is at 
the lower end, i.e. at the diffusion scale. A similar result was 
reported by Nakamura & Li (12007b who studied the effects of 
explosive partially collimated outflows on a background of de- 
caying turbulence and found a noticeable break in the velocity 
power spectrum at the outflow scale. Ossenkopf et al. ( 120011 ) 
focussed on self-gravity and demonstrated that self-gravity in- 
troduces a characteristic scale to the spectrum which evolves in 
time as more and more gas is accumulated in dense cores. This 
effect is probably only visible when using optically thin tracers 
sensitive for dense gas (Ossenkopf 12002b . which is not the case 
of our extinction maps. However, the observational data are not 
fully conclusive. While Ossenkopf et al. (12001b find a negative 
slope in Serpens (Testi et al. 11998b . a more recent study per- 
formed by Ossenkopf, Krips, Stutzki (2008b) based on a large 
mm-continuum map in pOph (Motte et al. ,1998j did not confirm 
this and instead reported the typical shape of the A-variance as 
discussed above. 
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Another recent application of the A- variance on model sim- 
ulations was performed by FedeiTath, Klessen, Schmidt (2009). 
They used the A-variance to characterise the density structure 
produced in supersonic, isothermal hydrodynamic models with 
two limiting cases of turbulence forcing: solenoidal (divergence- 
free) vs. compressive (curl-free) forcing. Table 2 shows that 
the compressive forcing models produce systematically higher 
values of fi than the solenoidal ones. All values are, however, 
at the upper limit of the range typically obtained from molec- 
ular cloud studies. Interestingly, a high value of fi is better 
consistent with the A-variance spectra determined from diffuse 
Galactic HI emission (see Table 1), and seems consistent with 
Kolmogorov scaling of incompressible turbulence (Falgarone et 
al. 2004). The transition to incompressible turbulence, however, 
is expected to occur close to the sonic scale at about 0.1 pc 
(Vazquez-Semadeni et al. 120031 Federrath et al. 1201 Oi l. Scales 
larger than that are clearly in the supersonic regime. The steep- 
ening of the power-law index /? towards large scales in HI clouds 
and in the numerical model with compressive forcing is thus un- 
likely to represent incompressible turbulence. It is more likely 
that the cloud-like structures on intermediate scales (see Fig. 1, 
right panel in Federrath et al. I2009I I produce a relative steepen- 
ing of the spectral index on large scales. A similar effect might 
occur for dense molecular cloud structures embedded in the dif- 
fuse atomic HI gas. In contrast, the steepening of fi on scales 
smaller than 0.1 pc (see Bensch et al. 2001) indicates a transi- 
tion to coherent cores (Goodman et al. 1998) with transonic to 
subsonic, almost incompressible turbulence, because this transi- 
tion is likely to occur close to the sonic scale (Vazquez-Semadeni 
et al. 120031 Federrath et al. l20T0l Pineda et al. |20T0] |. 

Table 2. Values of /? from hydrodynamic simulations of super- 
sonic, isothermal turbulence with solenodial (divergence-free) 
and compressive (curl-free) forcing (Federrath et al. l2009l l. 



Forcing p 


3d-hydrodynaimc simulation 


solenoidal 


2.89 


3d-hydrodynaimc simulation 


compressive 


3.44 


2d-projection 


solenoidal 


2.81 


2d-projection 


compressive 


3.37 



5. Summary 

We presented a 35 square degrees '^CO 1^0 molecular line sur- 
vey of Cygnus X, taken with the FCRAO, and visual extinction 
(Av) maps of 17 Galactic clouds, obtained from near-IR 2MASS 
data, in order to analyse the spatial structure of molecular clouds 
using the A-variance method. 

For Cygnus, we found no single characteristic scale in the 
Ay-map at all. A double-peak spectrum of the A-variance was 
revealed by '^CO 1— >0 data with peaks at 4 pc and 40 pc. The 
4 pc scale corresponds to the typical length of small filaments 
ubiquitously found in low-mass star-forming nearby clouds, 
while the second peak around 40 pc arises from the largest 
filamentary structures seen in Cygnus. This scale is below the 
whole extent of the cloud (~100 pc). An alternative explanation 
is that the 4 pc scale is the characteristic scale when the '^^CO 
1 — »0 line becomes optically thick (typically up to an extinction 
of about Av=5"'). Thus, '^CO observations may generally 
provide a good reflection of the structure measured only up 



to this Ay =5'" limit. The exact location of the peaks on the 
Delta-variance spectrum gives some measure for the turbulence 
in the molecular cloud. They indicate a transition from more 
coherent structures, where the '^CO emission falls into a narrow 
velocity interval, to a wide-spread turbulent distribution leading 
to a dilution of the optical depth in very broad lines. Applying 
the A-variance on velocity channel maps of molecular line data 
will help to understand better this point. The origin of energy 
injection in Cygnus can be manifold. It is known that the region 
is strongly affected by radiation from the Cygnus OB2/OB9 
clusters, expanding HII regions, a large number of outflow 
sources, and at least one supernova-remnant. 
The power spectral index fi for Cygnus , determined from the 
Delta-variance spectrum cr^ oc L^^^ with the size parameter L, 
has a value of j6=2.53 which is at the lower end of typical values 
for molecular clouds (typically between 2.5 and 3.3). 

The A-variance spectra obtained from the Ay maps show 
differences between low-mass star-forming (SF) clouds and 
massive giant molecular clouds (GMC) in terms of shape of 
the spectrum and values of the slope fi. Many of the low-mass 
SF clouds have a double-peak structure with characteristic size 
scales around 1 pc and 4 pc. scale at which a large-scale SN 
shock or expanding HII regions sweeping through the diffuse 
medium is broken at dense clouds, turning the systematic ve- 
locity into turbulence. GMCs show no characteristic scale in the 
Ay-maps which can be ascribed partly to a distance effect due 
to a larger line-of-sight (LOS) confusion. The values of fi show 
no clear trends, although there is a tendency that all fi values de- 
termined from '''CO are systematically higher than the ones ob- 
tained from the extinction maps and that high-mass SF regions 
have lower values of /? than low-mass SF regions. We do not see 
a trend with distance, despite of the fact that this might have been 
expected since shorter distances to a cloud resolve smaller scale 
structures and thus fi should be higher A comparison between 
the A-variance determined from model simulations (supersonic, 
isothermal hydrodynamic models with solenoidal (divergence- 
free) vs. compressive (curl-free) forcing by Federrath, Klessen, 
Schmidt ( I2009I I and observations shows that the model values 
are systematically higher, but still consistent with the observed 
range of A-variance slopes across different clouds. 

Overall, this study shows the full complexity of cloud struc- 
ture analysis, giving some indications of physics studied with 
the various methods. However, it leaves many open questions, 
which clearly require more work in comparing observations and 
turbulence models in the context of different methods. 
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Appendix A: Data qualtity 

The data points of the surveys were irregularly sampled and thus 
regridded to a fully sampled 22". 5 grid. The aliased noise power 
was minimized and the full resolution of the telescope retained. 
The convolution kernel used was 
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Fig. A.l. Rms map of the '^CO (C'^O) data in K. 



Znax cx/K,„ax 

with Ji the first-order Bessel function, x the distance from 
the observed data point, R,„„v the truncation radius {R,„ax-3), a 
= 0.9009, b = 0.21505, and c = 3.831706. n(/?„,„) is used as a 
pill box function with 11 = 1 for R<R,„ax and 11 = for R>R„,a t. 

The gridded spectra were converted and stored in CLASS 
spectra files. CLASS is a part of the GILDAS package and was 
used to further reduce the data. The data for small regions, typi- 
cally 30' X 30', were reduced in a data reduction pipeline. First, 
all spectra in the region were averaged to determine the win- 
dow around the emission to be excluded from the baseline fit. 
Second, a first to third order baseline was subtracted from all 
spectra. The reduced spectra were then written to 3-d data cubes 
(RA, DEC, velocity). Afterwards, cubes of the individual regions 
were merged to a single 3-d data cube for application of the A- 
variance. 

In order to quantify the quality of the survey, we produced 
rms noise maps of all line tracers, the one of '^CO is shown in 
Fig. lA.ll It becomes obvious that the noise is not uniform but 
shows a checkerboard structure and the edges are more noisy 
than the center. Both effects are due to the observing mode. In 
the OTP mode, a single position is passed by several pixels of 
the receiver array, resulting in a dense sampling and thus, an 
improved signal-to-noise ratio in the inner regions of an OTF 
map. Towards the map edges, the sampling is not as dense caus- 
ing higher rms values. Since the individual maps were observed 
with overlap to assemble the final large scale data set, the overlap 
regions are sampled a little more densely. The lower rms values 
in these overlaps create the checkerboard structure. The plot also 
highlights the areas which have been observed twice, namely the 
DR21 region. 



Appendix B: Extinction maps from dust reddening 
of 2I\/IASS sources: the example of Cygnus X 

The main purpose of developing AvMAP instead of using other 
programs (e.g. NICER; Lombardi et al. l2001l l is that it can han- 




Fig. B.l. Map of the extinction (Ay) toward the Cygnus X region 
obtained from the reddening of 2MASS point sources. The Ay 
derivation is optimized for a distance of 1.7 kpc, and is over- 
estimated for the more nearby complex of the North America 
and Pelican Nebulae. The Cygnus X complex is well visible and 
appears as mostly two large complexes referred to in Schneider 
et al. (I2OO6I 1 as Cygnus X North and South. One can also note 
that the S 106 cloud is well visible in the southern part of Cygnus 
X. 



die foreground stars up to larger distances and without assump- 
tions or off-field evaluation of the density of these foreground 
stars. Like that, we are able to image the most distant clouds. To 
illustrate our technique of deriving extinction maps from near- 
IR 2MASS data we describe in detail, as an example, the case of 
Cygnus X. 

The resulting Ay map for a field of 10° x 10° toward Cygnus 
X is displayed in Figure IbTI In this large field, a total of 4.5 x 10^ 
stars is found with an average of 12.4 stars per arcmin". The in- 
dividual Ay are measured as the quadratic, uncertainty weighted 
average of the two Ay estimates obtained from the two 2MASS 
colors [J-H] and [H-K]. The individual uncertainties are given 
by the quadratic sum of photometric uncertainties (from the cat- 
alog) and of the uncertainties on the intrinsic colors of a typical 
Galactic star (see for instance Lombardi et al. 2001 120011 for 
more details). The adopted values are [J-H]o - 0.45 + 0.15, and 
[H-K]o = 0.12+0.05 with the associated uncertainties measured 
from the color dispersions for a population of Galactic stars as 
measured using simulations with the Besangon Galactic models 
(Robin et al. l2U03l see below). 

B.1. Prediction of tlie density of foreground stars 

The above described method to estimate dust column densities 
entirely relies on the hypothesis that the near-lR sources are 
background stars. However, for a distant region like Cygnus X, 
a significant density of foreground stars is expected. To properly 
measure the dust column density, the foreground stars have to be 
filtered out (see detailed discussion in Cambresy et al. 120021) . 

We use the predictions of the Besan9on stellar models 
to derive a systematic, independent estimate of the expected 
number of foreground (nfg) stars for any Galactic direction. 
Catalogs of sources have been simulated using the online sim- 
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Fig. B.2. Variation of «fg measured from simulations using the Besangon model (see text) (a) in the Galactic plane as a function of 
longitude; and (b) in the Cygnus X region as a function of latitude. The variation with longitude is actually expressed as a function 
with rcc (see text). These predicted foreground star densities have been derived for a distance from Sun of 1.7 kpc. The adopted 
Unear interpolations are displayed as dashed lines. 



ulator (http :// model . obs-besancon.fr/) for a list of representative 
Galactic directions. For each simulated catalog, the number of 
stars brighter than the 2MASS completeness levels for all filters 
(J=17.0, H=16.5 and K=16.0), and at a distance smaller than the 
cloud distance (d = 1.7 kpc) is calculated to derive the expected 
density of foreground stars. Figure lBT2] displavs the resulting val- 
ues of nfg for 20 positions on the Galactic plane (lat=0) from 
lon=10 to 90°, and for 21 other positions at three fixed longi- 
tudes (lat=75°, 80 and 85°) from lat=0 to 10°. Note that nfg is 
expressed as a function of the Galactocentric distance rcc (down 
axis) using 7?^^ = 8.5 kpc. The Cygnus X region corresponds to 
rcc between 8.22 and 8.52 kpc where «fg is decreasing from 
~ 2 to 1.8 star per arcmin-, which is of the order of 10 % of 
the average density of 2MASS sources. The effect of latitude is 
larger with a typical decrease from ~ 1 .9 to 1.1 star per arcmin- 
for lat=10°. In order to filter out foreground stars, we then as- 
sumed that they are always the bluest stars in each direction 
which might not be always fully correct for the lowest extinc- 
tion regions. 



B.2. Distribution of young stars 

Young stars are usually more luminous than the population of 
background stars. If they are low or intermediate mass young 
stars, they are still above the main sequence, and if they are mas- 
sive, they are luminous and necessary young due to their fast 
stellar evolution. In the 2MASS catalog, these young stars are 
not necessary brighter than the surroundings due to the flux at- 
tenuation (extinction) of the parental cloud. However, follow- 
ing the same idea like to derive the average total extinction, the 
fluxes of all 2MASS sources can be de-reddened to evaluate their 
intrinsic brightness. Since the young stars form in clusters, maps 
of the distribution of the brightest IR sources after deredden- 
ing should trace well the embedded clusters. We thus use the 
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Fig. B.3. Map of the source density with a dereddened K band 
(2jum) magnitude smaller than 13. The linear scale is given as 
a slider on the right and is expressed in number of sources per 
arcmin^. The Cyg-OB2 association in the middle is clearly visi- 
ble (see also Knodlseder l2000l l. A number of embedded clusters 
inside Cygnus X North could be the seeds for a newly formed 
OB association. 



same Ay-mapping procedure to produce at the same time den- 
sity maps of the brightest dereddened stars. For the flux cut to be 
used, a good compromise between the 2MASS sensitivity and 
the typical brightness of young stars has been found to be be- 
tween K=12 and K=13. With a cut at K=12 a large number of 
cluster members can be found. In contrast, a cut at K= 1 3 is more 
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selective and can secure the detection of a cluster difficult to rec- 
ognize. 

Figure IB. 31 below displays the map of the density of stars 
with K>13 in the Cygnus X region after dereddening. It clearly 
shows a large-scale distribution of stars in the Galactic plane 
with several extended, higher concentration of stars. On the top 
of this large-scale distribution a number of small-scale clusters 
are apparent. A fraction of these clusters were also recognized by 
eye inspection of the 2MASS images by Dutra & Bica (1200 11 1, 
and using a simplified version of the present imaging procedure 
by LeDuigou & Knodlseder (l2002l l. 



B.3. "AvMAP", the fortran implementation of tine method 

In practice here for the Cygnus X field displayed in Fig. IB. H and 
IB. 31 the visual extinction Ay has been evaluated by averaging 
individual Ay for stars inside gaussian beams for each position 
in the map and after excluding the required number of proba- 
ble foreground sources. We adopted a grid of 450 x 450 pix- 
els with 80" spacing, and a resulting gaussian resolution of 1 .9' 
(Nyqvist sampling). The typical number of 2MASS sources per 
beam varies typically from 15 to 25 but goes down to close to 
10 in regions of high extinctions, and can even reach the typical 
number of foreground sources (4 to 5 per beam) for the darkest 
regions. This is particulary the case for the region at Ion ~ 84.7° 
which corresponds to the North America and Pelican Nebulae 
which are known to be at a distance of only 800 pc and for 
which, therefore, the number of foreground stars removed is too 
large leading to an over-estimate of Ay. 

In order to compute the different steps to derive Ay and the 
density of young stars in each pixel of large maps, we coded 
a FORTRAN90 program with an inner algorithm that is opti- 
mized to reduce the number of accesses to the (large) catalog 
of sources. In fact, for any large map, several millions of stars 
have to be processed with their coordinates, fluxes and uncer- 
tainties. It is impossible to go through millions of sources for 
each pixel when typically 10^ pixels are required to image the 
field. Instead, for each step of the calculation, each source is read 
only one time, and it is the contribution of this source to the final 
result which is calculated. Altogether, 3 paths/steps are required 
to compute Ay and the density of sources nK. 

The most recent version of AvMAP produces images directly 
in FITS format using the astrolib package. The best supported 
projection method is gnomonic even if some others implemented 
projection may work fine too (less tested). 

Appendix C: Figures of radio continuum and mid-IR 
emission of Cygnus X 

Figures C.l and C.2 show an overlay of '^CO 1— >0 emission 
from FCRAO on a radio continuum survey at 1420 MHz ob- 
tained from the CGPSH and on mid-lR emission at 8 jum from 
the MSX-satellite, respectively. The known thermal Hll-regions 
from the Downes & Rinehart survey ( |1966t are indicated. Since 
the mid-lR emission traces warm PAHs (Polycyclic Aromatic 
Hydrocarbonates) there is a very good correlation to the radio 
continuum of the Hll regions that serve as the heating sources 
for the dust. The molecular line emission is partly correlated 
with the mid-lR/radio continuum emission (bright-rimmed bor- 
ders of the cloud like for SI 06) but in other regions uncorre- 
cted (for example the emission seen at lower galactic latitudes). 
Examples for a stellar-wind/radiation driven bubble sre known 



sub— map 1 



aub— map 1 



Fig. D.l. Upper panels left: fBm-structure right: chess board. 
The three regions from which the A-variance (lower panels) was 
determined are indicated in the plot. 



http://www 1 .cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cgps 



for by NGC6913 (l=77°,b=0.5°)(also Schneider et al. l2007b and 
around the main OB cluster in Cygnus X (OB 2, l=80.2°,b=r). 
The nicely defined bubble at l=78.2°,b=2.1° is the well-known 
SNR gamma Cygni. 

Appendix D: Error calculation 

To test the significance of the A-variance, in particular the in- 
terpretation of the Poisson error bars on the interpretation of the 
spectra, we performed some numerical experiments to assess the 
actual impact of the "counting error" that is quantified in terms 
of the Poisson error bars. The counting error is due to the limited 
sampling of an overall (considered as infinite) structure by a fi- 
nite observation. We mimick this situation by selecting different 
parts of a larger periodic structure. 

Figure ID. II shows two different test maps, one fractal 
Brownian motion (fBm) structure and one periodic structure, a 
chess board. (See Ossenkopf, Krips, Stutzki ( I2008al l for more 
details.) The A-variance was determined from the whole map 
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Fig. C.l. The Cygnus X region in color scale at radio wavelengths (1420 MHz) from the Canadian Galactic Plane Survey. Contours 
outline the '^CO 1— >0 emission from the FCRAO. The stars indicate the most massive members of the Cyg OB2 association and 
crosses with numbers the HII regions DR4 to DR23. 



and from three regions with different sizes and is displayed in 
the lower panels. It becomes obvious that when observing differ- 
ent subregions of a larger structure, the A-variance spectrum of 
the substructure covers the full range of values indicated by the 
Poisson error bars, i.e. the application of the Poisson statistics 
for the computation of the uncertainty is confirmed. However, 
we find that the A-variance values at neighboring lags are not 
independent, but strongly correlated. The different view on the 
large structure "bends" the whole A-variance spectrum upwards 
or downwards, mainly at large lags, but it hardly introduces mu- 
tual variations between the A-variance values for similar lags. 
Consequently, we think that all the discussed structure in the 
A-variance spectra, including primary and secondary peaks, are 
significant because the error bars only describe the total uncer- 
tainty of the spectrum as a whole, but not variations of the Delta 
variance between neighboring lags. 
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